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^ ' Abstract. The systematic approach to study bound states in quantum chroniodynam- 

, ics is presented. The method utihzes nonperturbative flow equations in the confining 

^■f-^ ' background, that makes possible to perform perturbative renormahzation and to bring 

^SJ , the QCD Hamiltonian to a block-diagonal form with the number of quasiparticles con- 

serving in each block. The effective block-diagonal Hamiltonian provides constituent 
description for hadron observables. The renormalized to the second order effective 
Hamiltonian of gluodynamics in the Coulomb gauge is obtained at low energies. The 
\ masses for scalar and pseudoscalar glueballs are predicted. 
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§ : I INTRODUCTION 

r-| ■ One of the most difficult and less understood problems in quantum cliromody- 

Q-i! namics is the treatment of the bound state systems. There are different sources 

p of difficulties. For example, it is quite a common observation known in the spec- 

troscopy, that the splitting between the vector mesons does not depend on flavor, 
say 



m{p) — m{p) m{il)) — m{J /ijj) (1) 

H ; 

which is true experimentally. This fact can not be explained in terms of the canon- 
ical QCD interaction, which is given essentially by the strong coupling constant. 
Indeed, if expressed in terms of an invariant mass s, Eq.(l) implies that the J/ip 
is dual to a much larger interval of s than the p because the c-quark is heavy. 
However, the coupling constant runs as a function of an invariant mass, a(s), and 
is flavor blind, thus the canonical QCD interaction should be much weaker for the 
J/ip than for the p. This suggests that something is missing when described only 
in terms of the perturbation theory. The strong coupling constant alone does not 
provide for strong interactions being strong. 

Consider the scaling of Quantum Chromodynamics from high to low energies. In 
the ultraviolet region (at the bare cutoff scale A oo) the strong interactions are 
given by canonical QCD, which is conformally invariant, in particular this means 
scaling invariance (there is no scale in the theory) in the chiral limit. Moreover 



a perturbative treatment is possible due to asymptotic freedom. In asymptotic 
free theories (QCD) the couphng constant grows at low-energies and gets strong, 

that stops the asymptotic freedom at some moderate scale Aq ~ ^qcd- This scale 
appears in the theory when the perturbative rcnormalization of the coupling con- 
stant is performed, and the Aq is the Landau pole in the effective running coupling 
constant provided the rcnormalization group invariance. This is called dimensional 
transmutation, when scaling invariance breaks through the rcnormalization. Then 
experiment tells us that this is not the only scale in the theory. There are at least 
two more characteristic scales in the hadron physics. The mass gap of the hadron 
bound state, say, given by the square root of the string tension, where the nonper- 
turbative phenomenon of confinement takes place and the bound states of quarks 
and gluons form. The scalar vr-meson sets up a scale of chiral symmetry break- 
ing - another phenomenon of nonperturbative physics. The scales are displayed as 
m{n) -C < J^qcd- 

To summarize: Knowing perturbation theory alone is not enough to describe 
bound states in hadronic physics. 

Other difficulties in the QCD bound state problem are of a more general kind - 
relativistic nature of quantum field theory. First, the number of particles in any 
state is not fixed because of particle creation and annihilation in vacuum. Second, 
there are states with (infinitely) large energies. An attempt to treat both problems 
was done in [3] by using the method of flow equations. The idea is to find a unitary 
transformation that transforms the Hamiltonian operator to a block-diagonal form, 
where each block conserves the number of particles. The Hamiltonian matrix can 
be represented in the particle number space as 



where P and Q are projection operators on the subspaces with different particle 
number content. The flow equations [1] bring the Hamiltonian matrix Eq.(2) to 
the form 



where the two blocks of the effective Hamiltonian decouple from each other. It 
may be simpler then to solve for bound states within one block, say PHef/P, than 
to diagonalize the complete Hamiltonian, Eq.(2), of the original problem. Since, 
generally, the number of particles in P and Q spaces is arbitrary, one can reduce 
in this way the bound state problem with many particles to a few body problem. 

It turns out that the second question on possible presence of the ultraviolet diver- 
gences is solved also by the method of flow equations. Flow equations perform a set 
of (infinitesimal small) unitary transformations, where the flow parameter, which 
controls the transformation, has the dimension of the inverse of the energy square. 
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(3) 



I ~ 1/E'^. Therefore by using flow equations to block-diagonalize the Hamilto- 
nian one ehminates the particle number changing contributions not in one step but 
rather continuous for the different energy differences in sequence. This procedure 
enables one to separate the ultraviolet divergent contributions and to find the coun- 
terterms associated with these divergences. This covers the UV-renormalization for 
Hamiltonians [2]. 

To summarize: More generally, flow equations perform Hamiltonian renormal- 
ization in the "particle number" and in the "energy" spaces in the sense that the 
effects of the high Fock states and the effects of the large energies, respectively, are 
encoded in the effective low-energy Hamiltonian, which operates in the space of the 
few low Fock components. 



This program was applied quite successfully to QED to calculate the positronium 
spectrum [3]. The key is the validity of the perturbation theory in the bare coupling 
constant for the characteristic energy scale of positronium bound state. Obviously, 
it is a bad idea to apply naively the same scheme for QCD. It is not possible 
to find the fixed number representation for the Hamiltonian in the case of strong 
interactions, where one does not have any control over the process of the creation 
and the annihilation in vacuum of bare quarks and gluons with small current masses. 
In the language of flow equations convergence can not be achieved when calculated 
in terms of bare parameters. 

The way out is suggested by nature itself. One should consider "confined 
QCD". By using fiow equations one constructs then the effective QCD Hamil- 
tonian Hgff{q,g), where current quarks and gluons acquire masses of the order of 
TTiconstituent ~ IGeV and become constituent degrees of freedom (see below). The 
value of the constituent mass plays the role of the energy gap between the sectors 
in Hgff. Schematically, the block- diagonal effective Hamiltonian has the form 
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(4) 



where q and g are constituent quarks and gluons, respectively; empty cells denote 
zero, to the order calculations are done, matrix elements. To this order the different 
sectors of the effective Hamiltonian, Eq.(4), describe approximately, when going 
down, the bound states of mesons, hybrids, glueballs. Actually, such a description 
with the fixed number of constituents is quasiclassical and nonrelativistic, and is 
known from the constituent quark model. Physically, the picture is the following. 
The strong confining interaction, acting inside each "diagonal" (particle number 
conserving) sector, produces heavy gluon (quark) only in the small volume - in the 
"bag". No free propagating heavy gluons (quarks) are produced. Therefore the 
"bags" do not interact with each other, decoupling in i^e// ^-nd approximating the 
hadron bound states. 



The matrix elements of the "off-diagonal" (particle number changing) sectors 
are governed by the canonical interaction - typically by the Coulomb term of the 

strength equal to the inverse of Bohr radius or the current quark mass, rricurrent is 
of the order of several MeV. In the presence of the strong confining interaction in 
the "diagonal" sectors the mixing between the sectors is strongly suppressed. One 
can introduce a small parameter, say, 

Vi2 ^ TTT'current n t n nt /r:\ 

as— — ~ 1 • ~ 0.1 — 0.01, (5) 

El — E2 'ITT'constituent 

where V12 is the Coulomb interaction, and Ei — E2 is the energy difference between 
the first and the second "diagonal" sectors. Perturbation theory with respect to 
the small parameter, Eq.(5), holds between the sectors (but not inside the sector 
where the confining interaction is strong). By applying flow equations to block- 
diagonaUze the Hamiltonian one gets to leading order a closed chain of decoupled 
equations, which can be solved analytically. The whole is true provided there is a 
strong confining interaction in the "diagonal" sectors. 

To summarize: In the theory of strong interactions confinement is important to 
provide the bound states. In the present approach confinement makes it possible 
to bring, by fiow equations, the QCD Hamiltonian to a block-diagonal form with a 
fixed number of quasiparticles in each sector. The elementary degrees of freedom 
(quasiparticles) become constituent quarks and gluons, which acquire masses of or- 
der IGeV. The block-diagonal effective Hamiltonian approximately describes then 
the different hadronic bound states. 



The main idea of the approach is to find the representation for QCD Hamiltonian 
with the fixed number of quasiparticles, where the sectors with different particle 
number content decouple from each other. There can be some special cases when 
one should take into account the mixing between the sectors. In other words the 
physical state is not given by the pure component of the composite system. The 
mixing between the high excited state from the previous sector and the ground state 
from the next sector of the effective Hamiltonian may be possible (for example, the 
mixing between some excited meson and the low lying hybrid state). In systems 
with light quarks the influence of coupled channels can be essential. In the strongly 
coupled effective meson models one includes the effects from the coupled channels 
directly by mixing the scalar and pseudoscalar channels {qq and qqqq). The effect 
is about 50 percent. 

Special consideration is required in the case of the light quarks, where chiral 

symmetry breaking (CSB) is important. The present approach includes confine- 
ment and is like the "bag model" or the "constituent quark model", but it does 
not include CSB. By implementing CSB in this picture, the scalar 7r-meson can be 
viewed as a bound state of the two constituent quarks and simultaneously manifests 
the Goldstone nature. 

Motivation: In order to disentangle the both problems of confinement and CSB 
we consider the pure gluodynamics (see the next section) [5]. The motivation for 



this study is to set up a kind of a constituent gluon model, with the confining 
interaction imposed, to describe glueball bound states. 

The specific difficulty of QCD is that the canonical QCD Lagrangian docs not 
manifest explicitly confinement. As far as the mechanism of confinement is con- 
cerned one can proceed along several ways. To reveal confinement one uses the 
suitable formulation of QCD: lattice form, or the special choice of the gauge fix- 
ing (for example, maximal Abelian projection). Another option is to study other 
than QCD theories, but that have the same infrared behavior as QCD and are con- 
fined: Super Yang Mills theory, some toy gauge models (for example, Abelian Higgs 
Model). There may appear some unphysical degrees of freedom in these theories. 
If one is not interested in the mechanism of confinement, one includes the latter 
explicitly into QCD. The simplest way is to use the potential model, successfully 
tested in phenomenology, where the potential between the color charges is given by 
a sum of Coulomb and confining interactions. This suggest the definite choice of 
the gauge for the Hamiltonian. We work in the Coulomb gauge, where the Coulomb 
interaction arises from the gauge fixing procedure. We add then confinement to be 
able to block-diagonalize to the effective Hamiltonian, which describes the bound 
states. 

The Coulomb gauge is the natural gauge to get the constituent hadron picture. 
The Coulomb interaction appears there not as a perturbative (propagating in time) 
one gluon exchange, but rather as a solution of the gauge fixing constrains. There- 
fore the Coulomb term describes an instantaneous interaction, which is consistent 
with the nonpropagating massive gluon arising in our approach. Note, that the 
massive gluon mode arises only in the presence of confinement and confinement 
sets the scale for the gluon mass. The obvious drawback of the approach is the 
violation of the gauge invariance by the massive gauge fields. 

To summarize: The Coulomb gauge has an appealing property of the simple 
extension of the model to the confining case and is consistent with the constituent 
picture for hadrons. 

II LOW-ENERGY GLUODYNAMICS IN THE 

COULOMB GAUGE 

As noted above we applied the method to pure gluodynamics [5]. In this section 
we outline the strategy of this study. 

1. QCD Hcan- 

The starting point is the canonical QCD Hamiltonian (pure gluodynamics) in the 
Coulomb gauge V • A = 0: Hcani-^, n), where physical degrees of freedom are the 
transverse gauge fields A and their conjugate transverse momenta 11. 

2. H,,r,^H{g^0) + Oig) + Oig^). 

We expand the canonical QCD Hamiltonian in the Coulomb gauge perturbatively 
to the second order in the bare coupling constant. Then to the leading order 



the Faddeev- Popov determinant can be approximated by unity, that reduces the 
instantaneous term, arising from the gauge fixing, to the pure Coulomb interaction. 

3. Current (perturhaUve) basis and the trivial vacuum |0). 

We choose the trivial (pcrturbativc) vacuum |0) and construct the pcrturbative 
basis of free (current) particles: a'''(k)|0) creates one (pcrturbative) gluon with zero 
mass, i.e. the gluon energy u!]^ — |k|, etc., and the vacuum is defined as a\0) — 0. 

We express the canonical QCD Hamiltonian (section 2) in this pcrturbative Fock 
space, and normal order the result with respect to the trivial vacuum state |0). 
Denote the normal ordered canonical QCD Hamiltonian as -.Hcan'- 

4. Regularization and pcrturbative renormalization (scheme). 

The normal ordered Hamiltonian :Hcan' contains ultraviolet (UV) divergent terms 
(UV- divergent loop integrals). We regulate UV-divergences by the cutoff function 
f{q,A) (the explicit form of the cutoff function is specified further). This is the 
first time when we have introduced an energy scale in the theory - the bare cutoff 
A — > oo. To remove the cutoff sensitivity we renormalize the Hamiltonian by 
adding the counterterms associated with these divergences. Schematically, the 
renormalized Hamiltonian is written as Hreni^) = Hcan + ^Xcri-^), where SXcri-^) 
is a set of (unknown) counterterm operators, which we define further ^. 

5. Flow equations perturbatively. 

To find the explicit form of the counterterms and to scale down the Hamiltonian we 
run flow equations perturbatively. Also the form of the cutoff (regulating) function 
is specified by fiow equations. Generally, fiow equations define the prescription of 
regularization and make possible to perform the pcrturbative renormalization ^. 
Technically, since the Hamiltonian depends on the cutoff scale through the flow 
parameter, one flnds in the given order of perturbation theory (FT) the divergent 
part of the difference between the Hamiltonian operators given at the two scales, 
say (if(A2) — if(Ai)) with Aqcd -C A2 < Ai < A. One absorbs then these diver- 
gences in the counterterms - local operators with the symmetries of the canonical 
Hamiltonian, to provide the renormalization group invariance (called in the con- 
text of Hamiltonian renormalization "coupling coherence" [2]). This completes the 
procedure of renormalization, performed by flow equations, to this order. One can 
proceed in this way order by order in FT to find (all) the counterterms system- 
atically. Note, that it is enough to find the gradient of the Hamiltonian in the 
energy space to define the counterterms. Renormalization group invariance (RGI) 
insures, that the renormalized Hamiltonian preserves the form of the (original) 
canonical Hamiltonian, but only the coupling constants and the mass operators 
(that are usually classified as relevant and marginal operators in renormalization 
group sense) start to run with the cutoff scale. (We do not consider here, at least 

In the given (perturbative) basis this equation reads :-f/'ren(A); = '-Hcan- + :<5-^ct(A):, where 
" :" stands for normal ordering in the (perturbative) vacuum 

^) Flow equations perform a set of unitary transformations to block-diagonalize the Hamiltonian 
H{1, lo) = U~^{1, lo)H{lo)U{l, lo), where I is the flow parameter with the connection to the energy 
scale I = 1/A^, Iq is the initial value corresponding to the bare cutoff A introduced by the 
regularization before (section 4). 



in the few lowest orders of PT, possible irrelevant operators, that may cause new 
type of divergences than are carried by coupling constants and masses). 

Using flow equations we run the effective Hamiltonian downwards from the bare 
cutoff A to some intermediate scale Aq ~ ^qcd, where perturbation theory breaks 
down. Due to the RGI the "physical gluon" stays massless through this perturba- 
tive scaling. We can not proceed with flow equations perturbatively further. The 
result of this stage is the renormalized (to the second order of PT [5]), effective 
Hamiltonian, defined at some compositness scale Aq: :iJren(A, Aq): with Aq ~ ^qcd 
and bare cutoff A — > oo, and semicolon means normal-ordering in the trivial vac- 
uum |0). Though the renormalized Hamiltonian is obtained in the perturbative 
frame, it can be represented (regardless of the Fock basis) in terms of the fields A 
and n [5]. (It is a consequence of the RGI). We denote the resulting renormalized 
Hamiltonian at the scale Aq as Hreni-^, -^o) ■ 

6. Confinement. 

We introduce confinement as a linear rising potential, that enables to run flow 
equations "nonperturbatively" (see introduction) until complete diagonalization of 
the Hamiltonian. In the renormalization group sense this "spoils" the theory: there 
arises the massive gluon mode. But the presence of confinement is necessary to 
find the representation with a fixed number of quasiparticles (constituent massive 
gluons) for the effective Hamiltonian, which provides the constituent picture (see 
introduction). Confinement (string tension) sets the (hadron) scale for the gluon 
mass. 

The instantaneous interaction contains two pieces, the sum of the Coulomb and 
confining potentials. Denote the renormalized effective Hamiltonian with confine- 
ment embedded as iJe//(A,Ao). 

7. Constituent (nonperturbative) basis and the QCD vacuum 

As far as confinement is introduced the trivial vacuum |0) and the perturbative 
basis of free (current) particles, ur^ = |k|, define no longer the minimum ground 
state. Therefore, we introduce the (arbitrary) basis, where the gluon energy cu^ is 
kept unknown, and is defined further variationally. Correspondingly, the (nontriv- 
ial) QCD vacuum \Q) is defined as a\Q) = 0, and the Fock space of constituent 
particles is given: a^Q) creates the quasiparticle with the energy a;k, etc. ^. The 
renormalized effective Hamiltonian ife//(A, Aq) at the scale Aq (section 6), written 
through the physical fields A and 11 and having confinement, is decomposed in the 
trial (constituent) basis and normal-ordered with respect to QCD vacuum \Q). The 
unknown gluon energy is variational parameter in the calculations. We combine the 
terms in the effective Hamiltonian in each particle number sector according to the 

^) The change of basis from the (perturbative) current, cjk = |k|, to the (nonperturbative) con- 
stituent, with some Wk, can be written as Bogoluibov-Valatin (BV) transformation from the 
"old", a,a^, to the "new", a,a^ , operators: = ch(^kQ!k + sh(/>kQ:^j^ with BV angle 0k given 
by ch0k = ^/'2{\^k/uJk. + -^/i^k/k). The connection between the "old", |0), and the "new", 
vacuum states is given |f2) = exp J^k th0kttk«Lk) I* used in the work [4] to transform 
the QCD Hamiltonian into the constituent basis. 



power of coupling constant 0{g^) {n = 0, 1, 2) ^. In the absence of confinement the 
effective Hamiltonian preserves the form of canonical Hamiltonian due to RGI, with 

the proper change |k| — ^ cj^- In the presence of confinement the canonical form is 
violated by the second order terms in the effective Hamiltonian, which contribute 
higher orders 0{g^),etc. in flow equations. 

We aim to find the effective Hamiltonian after the scaling downwards from Aq to a 
hadron scale, say ^/a. Since the effective Hamiltonian preserves the canonical form 
at least to the second order, the "perturbative" terms obtained by flow equations in 
section 5 match the "nonperturbative" terms arising when applied flow equations 
to i/e//(A,Ao). 

We denote the effective Hamiltonian in constituent basis as ::Hgff{A, Aq)::, where 
stands for normal-ordering in the QCD vacuum. 

8. Flow equations in the confining background. 

We run flow equations in the conflning background to block-diagonalize the effective 
Hamiltonian ::ife//(A, Aq):: in the nonperturbative basis and to find consistently 
all the terms to the second order. Free Hamiltonian and conflning interaction are 
included in "diagonal" sector, the triple-gluon vertex forms "nondiagonal" sectors, 
that should be eliminated. We bring the Hamiltonian ::i7e//(A, Aq):: to a block- 
diagonal form, where diagonal blocks decouple from each other including the second 
order. The leading UV-behavior of the arising to the second order terms is cancelled 
by the mass counterterm. Generally, this approach allows to include perturbative 
QCD corrections into nonperturbative calculations of many-body techniques. The 
resulting block-diagonal effective, renormalized Hamiltonian is given at hadronic 
scale, ::ife//(A, y/a)::. For simphcity we denote it as i?e//- 

9. Gap equation (variational calculations). 

The requirement of block-diagonal form does not flx the effective Hamiltonian com- 
pletely. The remaining freedom to unitary transform inside of each block is flxed 
by minimizing the ground state (the vacuum expectation value of the effective 
Hamiltonian with respect to QCD vacuum) 

d{n\Heff\n)/dw = (6) 

to flnd the trial gluon energy uCk). The variational function c<j(k) is deflncd. As a 
result the gluon acquires a nonzero mass, m ~ 0.5GeV, and an iterative procedure 
of flow equations is performed with respect to the small parameter 1/m. The 
next following sector in the effective Hamiltonian is suppressed by this factor, that 
provides the convergence for the flow equations (see introduction). 
10. Solving for H^ff. 

There are two parameters in the method, the two scales: Aq, which defines the 
counterterms and regulates the perturbative radiative corrections to the effective 
Hamiltonian H^ff, and ^/a, where a is the string tension defining the nonpertur- 
bative confining potential. 



^) The higher order terms in the effective Hamiltonian are suppressed by the inverse powers of 
(heavy) gluon mass, which is of order of hadron scale (see introduction). 



Since the effective Hamiltonian is block-diagonal, one can solve for the bound 
states in any interesting sector (actually in the few lowest sectors). We solve for the 
glueball bound state in the two-body sector. The result is the glueball spectrum. 

Ill SUMMARY 

1. Renormalization 

Renormalization was performed to the second order. We combined the individ- 
ual counterterms in one- and zero-body sectors to the resulting mass countert- 
erm, written in the field representation (independent of the basis) 6Xct{^) = 
m^Tr / (ixA^(x) with = —^Nc^A"^ Remarkably, when the quark sector is 
added in the same fashion, the algebraic coefficient in the propagator correction 
reproduces the QCD /^-function. This particular feature of the Coulomb gauge 
supports our regularization prescription, which follows from flow equations. 

2. Glueball 

We specify the two parameters: the string tension is defined by the lattice calcula- 
tions a — 0.2GeV'^; the cutoff Aq is found from the gluon condensate to agree with 
the result of the sum rules, the condensate term is obtained (G^) ~ 1.3 • 10~'^GeV^ 
with Ao = 4:GeV. 

The solution of the gap equation can be parameterized as ci;(k) = k + 
m{0)exp{—k / k) , where the effective gluon mass is obtained m(0) = 0.90GeV and 
K = 0.95GeV. 

The glueball mass spectrum in scalar and pseudoscalar channels is given in 
Tamm-Dancoff approximation in [5] . Roughly, the mass of the lowest scalar glueball 
1760MeF, is twice of the effective gluon mass m(0). 
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